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ABSTRACT 

We report the discovery of multi-scale X-ray jets from the accreting neutron star X-ray binary, 
Circinus X-l. The bipolar outflows show wide opening angles and are spatially coincident with the 
radio jets seen in new high-resolution radio images of the region. The morphology of the emission 
regions suggests that the jets from Circinus X-l are running into a terminal shock with the interstellar 
medium, as is seen in powerful radio galaxies. This and other observations indicate that the jets have 
a wide opening angle, suggesting that the jets are either not very well collimated or precessing. We 
interpret the spectra from the shocks as cooled synchrotron emission and derive a cooling age of 
~ 1600 yr. This allows us to constrain the jet power to be 3 x 10 35 ergs _1 < Pj ct < 2 x 10 37 ergs _1 , 
making this one of a few microquasars with a direct measurement of its jet power and the only known 
microquasar that exhibits stationary large-scale X-ray emission. 

Subject headings: ISM: jets and outflows, X-rays: binaries, X-rays: individual (Circinus X-l) 



1. INTRODUCTION 

Circinus X-l is an unusual X-ray binary (XRB). 
I Stewart et al.l (|1993f ) first noticed that Circinus X-l re- 
sides inside of a parsec-scale radio nebula inflated by 
curved jets. More-recent radio observati ons have ob- 
serve d jet emission on multiple scales ()Tudose et al.l 
120081 and references therein). Additionally, at the peak 
of the 40+ year X- r ay light curv e (see Figure 1 of 
iParkinson et al.ll2003| ). [Fender et all (|2004h detected ev- 
idence of superluminal motion from ultra-relativistic jets 
within a few arcseconds of the XRB, which would make 
this the only known accreting neutron star with an ultra- 
relativistic jet. With its resolved, parsec-scale radio jets, 
Circinus X-l has become an important stepping stone 
in our understanding of microquasars over the past few 
decades. 

Circinus X-l has been extensively studied at X-ray 
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w ayelengths. Type-I X -ray bursts, origi nally observed 
bvlTennant et al. (1986) and seen again bv lLinares et al.l 
(2010) during a recent flare, now firmly establish that the 
compact object is an accreting neutron star. Multiple 
Chandra gratings observations of the point source have 
revealed variable X-ray P-Cygni line profiles that were 
interpreted as high- velocity ou tflows (iBrandt fc Schuld 
l200a ISchuIz fc Brandtj[2002T) . Ilaria et all (120081) have 
claimed possible detection o f a highly-inclined precess- 
ing jet (similar to SS433, e.g.. lLopez et al.ll2006l ) through 
Dopplcr-shiftcd X-ray line emission. However, this in- 
terpretation is not unique as the emission is also con- 
sistent wi th simple orbital m otion of the neutron star. 
Moreover. ISchulz et alj (|2008D also analyzed that partic- 
ular data set and, with systematic uncertainties included, 
find no shifts. Finally, in the longest gratings observa- 
tion of the point so urce (zero-order, 50 ks) taken in 2005, 
iHeinz et"aT] (|2007f l discovered faint, diffuse X-ray emis- 
sion roughly coincident with the arcminute-scale radio 
jets. However, the observation had insufficient signal- 
to-noise to constrain the nature of the emission spec- 
troscopically. Similar diffuse emissio n was subsequentl y 
seen in a 50 ks HRC-I observation (Sol eri et al.1 12009). 
but no spectral information was available with this in- 
strument. 

In this Letter, we report on a follow-up deep Chandra 
im aging observ a tion o f the diffuse X-ray emission found 
by IHeinz et alj (|2007[ ) . We outline our observations in 
Section 2 and present an initial analysis of the emission 
in Section 3. Section 4 discusses the physical implica- 
tions. Finally, we summarize our results in Section 5. 
Throughout this Le tter, we assume a distance of 7.8 kpc 
(| Jonker et al.ll2007t) . which is in the middle of a wide 
range of distance estimates (4.1-11.8 kpc; Iaria et al. 
12503) . 

2. OBSERVATIONS 

We observed Circinus X-l on the S3 chip of the Ad- 
vanced CCD Imaging Spectrometer (ACIS) aboard the 
Chandra X-ray Observatory on 2009 May 1, in a contin- 
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Fig. 1. — Reduced event file showing the point source and the 
bipolar "caps" created by the outflow. Overlaid are contours from 
the lower resolution radio image (logarithmically spaced from 1.9 
to 25 mjy beam ; beam size: 18."1 X 16."5) from lTudose et al.l 
(2006). 

uous 99 ks exposure. Data were taken in timed exposure 
mode and telemetered in Faint mode. Data reduction 
and analysis were completed using CIAO and Shcrpa 
versions 4.2, XS pec 12.5.1, a nd ACIS Extract version 
2010-02-26 fAE: iBroos eFail [2010h . Compared to pre- 
vious gratings and HRC-I observations, this observation 
proved about an order of magnitude more sensitive to 
diffuse flux not only because of the long exposure and 
the lack of gratings, but also because the contaminating 
emission from the point source was at an exceptionally 
low level (Section 3.1). 

Radio comparison data were constructed from two 
separate observa tions. Standard calibrations using the 
Miriad software (|Sault et al.lll995T ) were applied to both 
sets of observations. A low angular resolution radio im- 
age was derived from data taken on 2001 August 3 dur- 
ing an 11 hr run at 1.4 GHz with the Australia Tele- 
scope Compact Array (ATCA) in 1.5A array configura- 
tion. Thi s observation a nd th e resulting images were pre- 
sented in lTudose et al.l (|2006l ) and are shown as contours 
in Figure [T] 

A higher angular resolution radio image (Figure [2]) 
was created from ~ 80 hr of observations spread from 
2009 December 30 to 2010 January 8, using the ATCA- 
Compact Array Broadband Backend in 6A configuration 
at 5.5 GHz (D. E. Calvelo et al., in preparation). 

3. ANALYSIS 

The reduced event file of the X-ray observation is 
shown in Figure Q] with a contour overlay of the low- 
resolution radio image of the large-scale radio nebula of 
Circinus X-l. This image clearly shows the existence of 
two bright emission regions ~ 30" from the point source. 

In order to highlight the morphology of the diffuse 
emission closer to the point source and to allow detailed 
comparison with the high-resolution radio image, we cre- 
ated an exposure-corrected, background-subtracted im- 
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Fig. 2. — Exposure-corrected and azimuthally-smoothed, 
background-subtracted X-ray image. Overlaid are contours from 
the higher resolution radio image (levels: -1, 2, 4, 8, 16, 32 times 
rms noise of ~ 6 /ijy beam -1 ; beam size: 3."02 X 2."21; D. E. 
Calvelo et al., in preparation). The point source, inner emission, 
and the X-ray caps described in Sections 3.1-3.3 are all evident 
here and well-matched to the radio contours. 

age (Figured]). Background surface brightness profiles 
were constructed by azimuthally averaging over nested 
partial annuli in the NE and SW quadrants (visually 
chosen to exclude the excess diffuse emission) and then 
subtracted from the exposure-corrected image. 

Finally, we note that on the largest scales, we find 
roughly axisymmetric soft emission from what we in- 
terpret as the dust-scattering halo, centered on the bi- 
nary, best visible in the adaptively-smoothed image in 
Figure El An in-depth discussion of the dust-scattering 
halo will be presented in a follow-up paper. 

Image analysis and X-ray-radio comparisons are pre- 
sented below, ordered from smaller to larger scales. All 
well-determined background/foreground point sources 
were masked and excluded from our analysis. 

3.1. The Point Source 

We caught Circinus X-l at the lowest observed flux to 
date (-Fn.5-8keV — 1 x 10~ n erg cm -2 s -1 ). Even so, the 
source is still significantly piled up with a pileup fraction 
of ~ 40%. We will present detailed model fits to the 
point source spectrum in a follow-up paper. 

3.2. An Arcsecond Jet 

Because the point source was so faint, the point spread 
function (PSF) from the XRB only dominates the in- 
ner ~ 5". This enabled us to search for jet emission on 
scales of a fraction o f a parsec, much sm alle r than what 
had b een possible in lHeinz et al.l (|2007| ) and lSoleri et al.l 
( 2009]). We found two sources within ~ 15" of the point 
source, which are both evident in Figure [2] 

1. In the W-NW direction (between position angles 
0° and +45°, measured counterclockwise from due 
W), we find a clear surface brightness enhancement 
separated from the XRB by ~ 5" and extending out 
to ~ 8". The emission appears to be resolved and 



3 




Fig. 3. — Adaptively smoothed three-color image ~ 4.'5 X 4/5 
(red: 0.5-2.9 keV, green: 2.9-4.2 keV, and blue: 4.2-7.0 keV). The 
large-scale dust-scattering halo from the point source (red) and the 
hard wings of the PSF (blue) are clearly distinct in color from the 
cap emission. North is up and east is left. 

is coincident with a source in the high-resolution 
radio image. No IR counterpart is found in ei- 
ther 24 fim 

MIPSGAI0 or 3.6-8/Ltm GLIMPSE^ 
Spitzer images. Given that the source is extended 
and aligned with the large-scale radio jet and with 
the X-ray caps reported below, we interpret this 
source as genuine jet emission. 

2. In the SE direction (at position angle 230°), we 
detect a bright, unresolved point source at ~ 13" 
from the XRB. This emission is coincident with a 
point-like radio source. It is spatially coincident 
with a very bright IR point source, with no opti- 
cal counterpart detected in an archival 5 minute 
HST/WFPC I exposure. 

From the FIRST (Faint Images of the Radio Sky 
at Twenty-Centimeters) log TV-log S, we estimate 
that the probability of finding a background radio 
point source of the same or larger flux within 13 
arcseconds from Circinus X-l is ~ 1%. In addi- 
tion, given the IR source density in the region, the 
likelihood of finding an IR background source with 
equal or larger flux coincident with the radio and 
X-ray point source within the MIPS PSF is also of 
order 1%. Neither scenario can be ruled out sta- 
tistically. Thus, we cannot draw any conclusions 
about the nature of the SE point source without 
further observational evidence, e.g., from optical 
spectroscopy, variability, or proper motion. 

A detailed PSF subtraction of the binary emission 
(complicated by effects of pileup) and spectral fits of 

13 http:/ /mipsgal. ipac.caltech.edu/ 

14 http:/ /www. astro.wisc.edu/sirtf/ 



these sources are beyond the scope of this Letter and 
will be presented in a follow-up paper. 

3.3. Extended Diffuse Emission: X-ray Caps 

The most obvious features in Figures QJ3] are the two 
diffuse emission regions between ~ 20" and 50" from 
Circinus X-l in the NW and SE directions. The po- 
sitions of both regions are con sistent with the d iffuse 
emission tentative ly reported in iHeinz et al.l (|2007f ) and 
iSoleri et"aLl ((2009). Based on their morphological ap- 
pearance (NW cap appears slightly concave) and their 
placement away from the binary along the jet axis, we 
will refer to both regions as "caps" throughout the rest 
of this Letter. However, the geometry of these regions 
could be complicated by projection effects, and other in- 
terpretations are possible. 

The caps have distinctly different X-ray colors from 
both the much redder large-scale diffuse background 
(dominated by the dust-scattering halo) and the dis- 
tinctly bluer PSF, as can be seen from the color image in 
Figure [3] The two caps have very similar surface bright- 
nesses (~ 2 x 10~ 16 erg cm~ 2 s _1 arcsec~ 2 ) and appear 
to be similar in their angular extent (NW: — 5°-65°; SE: 
185°-245°). However, they are asymmetric in their ra- 
dial extent (NW: - 22"-50", SE: ~18-45"). The inner 
NW jet emission (Section 3.2) has similar position and 
opening angles to the NW cap. 

The correspondence between the high-resolution radio 
contours (> 2a) and the X-ray cap emission in Figure[2]is 
striking. While the X-ray image has significantly higher 
angular resolution than the radio, it is clear that the X- 
ray caps align closely with the extended jet emission. 

The diffuse emission to the SE of the XRB in the radio 
image shows a sharp surface brightness drop at the posi- 
tion of the cap and what has been interpreted as a bend 
in the jet direction at this position, trailing to the south 
(referred to as "knot B" in lTudose et~aT]|2006[ ). We see 
no clear corresponding surface brightness enhancement 
to this trail in the X-rays. However, our sensitivity is 
compromised by the fact that the read streak runs right 
along this feature. 

We used partial annuli regions to extract the cap spec- 
tra. Because the background is dominated by the dust- 
scattering halo, the surface brightness of which depends 
on the angular distance to Circinus X-l, we carefully 
selected background partial annuli at the same radial ex- 
tent as the source regions. We recover a total of 6240 
source counts over 2860 expected background counts in 
both cap regions combined (0.5-9.5 keV). The resulting 
spectra are shown in Figure 01 

Initial extractions of each of the caps indicated that 
the spectral parameters are statistically indistinguish- 
able. Therefore, we jointly fit the spectra. We restricted 
our spectral analysis to 0.5-9.5 keV. Quoted uncertain- 
ties correspond to 90% confidence intervals. The spectra 
are well fit by an absorbed powerlaw with a neutral hy- 
drogen column of TVh = 2.24^28 x 10 22 cm~ 2 and a pho- 
ton index, V = 1.99±£|i (% 2 = 243.1 over 252 degrees 
of freedom), giving 0.5-9.5 keV cap fluxes of 3.2 x 10~ 13 
erg cm" 2 s" 1 (NW) and 2.7 x 10~ 13 erg cm -2 s" 1 (SE). 

The caps can also be fit by a thermal (APEC) model 
(x 2 = 237.7 over 251 degrees of freedom). The required 
absorbing column is very similar to the power-law case 
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Fig. 4. — Top: spectrum of the two large-scale caps and best-fit 
absorbed power-law fit (red: NW cap; green: SE cap). Bottom: 
deviations from the best-fit model. 

with N H = L93±o'.i8 x 1q22 cm" 2 . This is consistent with 
the fact that there is no obvious line emission present in 
the spectra, the APEC fit requires a low abundance of 
Z = OAltoir and a hi g h temperature of kT = 6.6t 2 ;g 
keV. When forced to solar abundance, the fit requires 
even higher temperatures. The 0.5-9.5 keV cap fluxes 
are 3.2 x KT 13 erg cm" 2 s" 1 (NW) and 2.7 x 1CT 13 erg 
cm" 2 s- 1 (SE). 

4. DISCUSSION 

The limb-brightened morphology of the X-ray image 
(which has significantly higher angular resolution than 
either radio image) suggests that the emission from the 
two caps originates at a shock in the outflow from Circi- 
nus X-l, seen in projection. This suggestion is supported 
by a sharp drop in surface brightness in radio and X-ray 
just outside of the caps and by the fact that the SE ra- 
dio jet seems to be changing direction at the location of 
the X-ray cap. The fact that the caps appear inside the 
large-scale radio nebula would be a result of foreshort- 
ening, since the jet axis is likely inclined with respect 
to the line of sight (though the actual inclination is un- 
known). The outflow itself appears largely X-ray dark 
(except for the inner NW diffuse radio-X-ray feature), 
similar to the X-ray cavities observed in many clusters 
with central radio galaxies. 

The caps span projected half-opening angles of 35° 
(NW) and 30° (SE). This implies that either: 

1. the outflow is, in fact, much wider than what would 
typically be considered a jet and might thus be bet- 
ter characterized as a non-thermal wind, or 

2. the jets are precessing with a fair ly wide open - 
ing angle, as is the case for SS433 (|Margonl fl984). 
and/or 

3. the jet or precession cone axis is very close to the 
line of sight, causing significant foreshortening of 
an intrinsically narrow jet. 

The wide opening angle of the inner NW jet would be 
consistent with (1) and (3), but to be consistent with (2), 
it would require a precession period of the jet that is short 



compared to the travel time through the X-ray emission 
region to blend the jet emission (as seen in the a rcsecond 
radio jets of SS433 lHiellming fc Johnston! 1 19811 ) . 

The claim of a h ighly relativistic jet (rj Gt ~ 16, 6 < 5°; 
iFender et alJl200"4] ) coupled with these new observations 
would require that such a jet has a precession cone 
that is within 5° of the line-of-sight such that the jet 
sometimes points very close to the line-of-sight. Such 
a geometry would imply a physical scale of ~ 20 pc 
from cap to cap. For comparison, this is an order of 
magnitude larger than the project ed distance b e tween 
the hotspots of XTE J1550-564 dCorbel et alJ l200"l 
and H1743-322 ([Corbel et al.ll2005l) but would still fit 
comfortably within the radio nebula surrounding SS433 
(jDubner et al.|[l998h . 

Below, we will discuss two possible interpretations of 
the radiative origin of the cap emission. 

4.1. Cooled Synchrotron Model 

The most likely interpretation of the spectra is syn- 
chrotron emission. In this case, the emission arises 
from the reaccelerated jet particles entering the termi- 
nal shock, as is obse rved in the hot spots of m any FR II 
radio galaxies (e.g.. iMeisenheimer et a.1.1 11989T ) . The X- 
ray emission in this case should be co-spatial with the 
radio, which appears consistent with the radio contours 
in Figure [2] 

The radio synchrotron spectrum of the large-scale neb- 
ula has a typ ical spectral index of or = 0.53 (with 
F v cx is~ aR ) (jTudose et al.ll2006f ). consistent with the 
standard power-law slope of first-order Fermi accelera- 
tion. For the canonical model of continuously injected 
power-law electrons one would expect a spectral index 
of ax = cur + 0.53 = 1.03 above the cooling break at 
frequency Vb, perfectly consistent with the observed X- 
ray spectral index of ax = 0.99. From extrapolation of 
the radio and X-ray spectra, we estimate the break fre- 
quency from uncooled to cooled synchrotron emission to 
be v\> ~ 2 x 10 16 Hz, which is uncertain by about an 
order of magnitude when uncertainties in the radio and 
X-ray spectral indices are used primarily because the X- 
ray and radio power-law slopes are extrapolated over a 
very large range of frequencies. 

Taking the emission regions to be spherical to low- 
est order, we estimate the equipartition magnetic field 
strength to be B cq ~ 50 /iG, which gives a minimum to- 
tal internal energy for each emission region of ~ 8 x 10 45 
erg. From the estimate of the break frequency, we esti- 
mate a cooling time, r coo i ~ 1600 yr. A different emission 
geometry will change these numbers by factors of order 
unity. 

Together, these numbers provide a robust lower limit 
on the total jet power of P c - lr > 3 x 10 35 ergs -1 , needed 
just to put the synchrotron emitting particles in place 
(not including any plasma currently in the jet or the en- 
ergy needed to inflate the large-scale radio nebula). De- 
parting from equipartition can only increase this number. 
This limit is independent of the assumed geometry and 
size of the emission regions. 

Conversely, we can place a rough upper limit on the 
outflow power from the fact that the large-scale radio 
nebula must be at least as old as T CO o i. U sing the same 
argum ents as presented in[Heinz (2002]) and lTudose et al.l 
( 2006) for the expansion of the radio nebula we find that 
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-Pcir < 2 x 10 37 ergs _1 , though this limit is not as robust 
since the magnetic field could well be out of equiparti- 
tion and the nebula could be significantly elongated along 
the line-of-sight, both of which would increase the total 
power. 

These limits on P clI confirm that the accreting com- 
pact object in this system is driving powerful jets into 
the interstellar medium (ISM). Taking the high 42 year 
average luminosity of Circinus X-l at about 80% of 
the Eddington lumino sity for a 1.4 Mq neutron star 
(jParkinson et al.l I2003D at face value would imply that 
the long-term average jet power is between 2% and 10% 
of the current average radiative power. 

Given that the average long-term accretion rate cannot 
greatly exceed the Eddington rate, we can derive a robust 
lower limit on the efficiency with which accretion power 
is converted into outflow power of T7j Gt = P C h /rnc 2 > 
0.034%^u. (see also iHeinz et al.l [20071 and ISoleri etafl 
[2009h . 



4.2. Thermal Model 

Given the equally good fit we achieved with a ther- 
mal emission model, it is worth discussing the alterna- 
tive scenario that the X-ray caps are the shocked ISM 
pushed ahead of the expanding radio outflow (similar to 
the shells around many X-ray cavities in gal axy clusters 
and around the microquasar Cygnus X-l; IGallo et al.l 
l200Tih . 

In the thermal model the X-rays should be located fur- 
ther out from the binary than the radio emission. How- 
ever, the high-resolution radio data indicate that the ra- 
dio and the X-ray emission are co-spatial. This requires 
a very low inclination angle in order to align the two 
regions by projection (which is consistent with the incli- 
nation required by the ultra-relativistic flow claimed in 
iFender et ai][200l . 

From the best-fit temperature of 6.6 keV, we directly 



infer a shock velocity of ~ 2.4 x 10 3 km s _1 . From the 
projected distance, we can infer an estimated travel time 
of ~ 500 / sin(0) years out to the observed shock positions 
for each of the caps. Again approximating the emission 
region as spherical, we estimate the mass and the total 
(thermal plus kinetic) energy of the shocked material for 
each of the caps to be ~ 2 x 10 32 g and ~ 2 x 10 49 erg, 
respectively. Combined with the shock travel age, this 
gives a minimum outflow power for both caps combined 
of P cir > 6 x 10 38 x sin (6) erg -1 s -1 , just to supply 
the X-ray emitting material (not including the putative 
power required to inflate the large-scale radio nebula). 
Unless the inclination angle is very small, this power is 
exceedingly large. Even for inclinations as low as those 
inferre d from the relativistic jet claimed in IFender et al.l 
( 2004) , the minimum average power would still be P c ; r > 
5 x 10 37 ergs _1 . 

While we cannot completely rule out a thermal ori- 
gin, the very high power required by the thermal model 
combined with the perfect spectral correspondence to a 
classic broken synchrotron spectrum and the spatial cor- 
respondence with the radio suggest that the synchrotron 
model is the more natural choice. 

5. SUMMARY 

We have presented an initial analysis of a deep Chan- 
dra imaging observation of Circinus X-l. We detect two 
diffuse X-ray caps that are likely the terminal shocks of 
powerful jets running into the ISM. In addition, we find 
an arcsecond-scale outflow between the XRB and one of 
the X-ray caps, coincident with the radio jet. These dis- 
coveries make Circinus X-l the first microquasar with 
both an X-ray jet and two stationary X-ray shocks, and 
one of only a handful of microquasars with a direct esti- 
mate of the jet power. 

S.H. and P.S. acknowledge support from NASA grant 
GO9-0056X. 
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